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1.  Introduction 


Covariance  and  correlation  natricea  and  their  principal  eoupcnsnts 
arc  basic  alaMota  in  noltivarlata  data  analysis.  Vat,  tbs  classical 
estiantes  of  thasa  quantities  ara  vary  s await! as  to  outliers.  In  i so sat 
yaara,  aavaral  robust  altasaatlvaa  bars  boon  propuasd  and  stud  lad  (see 
111,  (II,  |4|,  Ml,  (91,  HOI.  |14|.  |1S|.  and  |1*|).  Tbs  protoleo  aaM 
to  ba  that  no  om  ia  rnally  satisfactory.  Ha  natrix  elenant  Mtbsds  ara 
not  orthogonally  aqoirariaat,  da  net  mm  ywrwisa  to  prodaon  a  poaitiva 
dafinita  natrix,  and  ara  vary  difficult  to  study  theoretically,  Mila 
af finely  squivariant  R-eetiaatea  of  conrianoa  natrix  bars  quite  poor 
bresMnau  properties  in  biqh  dlnsnalon. 

Ibis  pap*',  together  with  (21 ,  proposes  and  diecuasaa  a  n aw  type  of 
aatiautor  for  cover  ience/oorrelat  ion  an  trices  and  tbeir  principal  esapon- 
Mta.  It  uses  a  Projection  Pursuit  IPP)  procotote  Isos  ID,  |5|,  Ml,  111, 
(HI,  1 1 2 1 >  with  a  robust  satinets  of  seal#  as  tbs  projection  index.  So 
we  call  it  the  ROBUST  rP-BSTIWMUR.  Its  idea  wax  originally  raised  by 
Rater  (sea  (lot,  p.  200  and  pp.  203-204). 

Tbs  pp  procedure  duals  with  hiqh  dtnanxlonal  dstai  it  aaarebaa  iow- 
djas ns jewel  projections  ubich  aaxiniza  (win ini as)  an  objoctivo  function 
eallod  PROJECTION  HBO. 
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m  followai  let  te  a  voteat  astieator  for  scale,  which  la  nasally 
waafcly  eoetianona,  and  a  €  B*  ten  vector.  Denote  the  diatrltetioa 
fsnetian  of  «f«  by  »?,  or  f(aVI)  vteo  needed.  The  eat  law  tore  for  prlaci- 
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DUm  that  Ima  M  cm  an  mwmm  IWI  Uni*  of  W  ainea 

«•  Ml—  tlMt  »•)  ia  aaakly  oMUama,  Uw  for  aar  F,  »(F?)  la  a 
c—laaoao  function  of  a  (a—  La— a  4.4  U>  Mlw),  —  tka—  alwaya  aalat 
at  (1  <  i  <  p)  «Meh  raaofc  tka  —I—  tnlal—n |  nl—  la  11.1)  1(1.1*11 
—  tka  oorroapcnftwy  raniona,  gaapaetloaly. 
okrlaaaiy. 
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Niwn  4^  la  a  probability  aoaaare  putting  a  polataaao  at  a. 
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it  is  aaay  to  show  that  If  Ixl  <  «,  thw 
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Oawamlr,  If  c  >  ,  than 
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On  tha  othar  hand, 
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Aa  on  will  aaa  In  tha  and  of  Sactloa  9.  it  la  lapartiat  that  S(r)  doaa  aof 
•labahava  aa  a  projection  indon  for  principal  oooponant  aatlnation.  that 
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1.  TilnrUBct 


Tha  study  of  robust  astlnators  of  oovuUbm  powa  ■  problani  Rw  to 
protUi  «  covar  lanes  salinate*  sblcb  hahavaa  utter  eoordlasts  ahanpss  ttes 
hh  usy  as  tha  classical  sstlaator  doss.  It  is  important  to  sstateliate 
sgslvar  lanes  propart las  of  a  nw  sstlaator  in  octet  to  tew  that  it  is, 
truly,  an  sat  lam  tor  of  dispart  ion.  Hilo  auction  about  that  tha  robust  PP- 
satinstora  do  Hava  squisarlanos  pcopartiss. 

ta  It  la  usntlonad  in  Suction  I  that  this  paper  oonsidsrs  tbs  purs 
dlapsralon  problsn,  i.a..  aasns  that  location  is  hnoan  and  f inad  at  0. 
Also,  wa  aaouno  throughout  that  »(*)  is  wakly  continuous. 

Mllbi  iht  »,(.)  al|nM  <k1w.1i  MwIkI.  n  mm  tkw 
to  denote  any  version  of  these  quantities. 


iwteOMH  3.1.  hobust  PS-sat  inatas  J^fri  (1  -  1,  2,  ....  pi 

and  $tr)  are  all  orthogonal  ago! variant.  I.a..  lot  0  ha  any  urtteo  panel 
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an  indeed  orthogonal  eqoi  variant,  jeet  lilta  ola— leal  apt  teat  oca. 

Nvm*  that  ru(?)  balance  to  a  p  dl-glen  elliptic  pm  dtillty 
daaaity  f  rally,  I. a..  FU.Vl  haa  a  density  fU/n 
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4.  Consistency  at  Elliptic  Probability  Penalty  fanlly 

g—  the  data  com  frxm  an  elliptic  probability  diatr  Unit  ion,  for 
which  pwtewTianet  haa  an  In  tall  lira  Interpretation  -  the  shape  of  the 
wwWr  lying  allipaa,  It  la  poaalbla  to  atm*  that  the  rotoest  PP-eet instore 
fin  conalatant  estinates.  Ihe  Idea  of  the  proof  la  alnplei  It  ooafeinea 
the  continuity  property  of  the  projection  late  l(')  (dlecneeed  in 
Section  3  tor  M-eatlnataa)  with  a  coapectneea  arqonent .  tat  putting  this 
into  operation  In  high  dinenaion  and  for  any  enderlylng  paeudo  covariance 
ia  quite  ccnpllcated.  To  cope  with  this  probien,  in  introduce  momm  lanes 
first,  then  start  with  sons  special  underlying  paando  covariances  to  achieve 
general  consistency  results. 
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This  paper  proposes  and  discusses  the  ROBUST  PROJECTION  PURSUIT  ESTIMATOR 
for  dispersion  matrices  and  their  principal  components.  This  estimator  finds 
robust  principal  components  by  searching,  successively,  for  directions  which 
maximize  (minimize)  a  robust  estimate  of  scale;  the  estimate  of  the  dispersion 
matrix  is  constructed  from  the  estimated  principal  components. 

These  estimators  are  shown  below  (under  mild  conditions)  to  have  a  number  o 
desirable  properties.  They  are  orthogonally  equivariant  and,  within  any 
elliptic  underlying  density  family,  asymptotically  affinely  equivariant. 
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Furthermore,  at  elliptic  densities,  they  are  consistent  and  weakly 
continuous  (l.e.,  qualitatively  robust).  Finally  they  have  good  quantitative 
robustness— their  breakdown  point  can  be  as  high  as  1/2. 

The  robust  projection  pursuit  approach  is  a  promising  alternative  to 
other  estimators  of  dispersion  matrices. 
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